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ABSTRACT 

in 

O 

| The prompt localization of gamma-ray burst (GRB) 050525A by Swift allowed 

the rapid follow-up of the afterglow. The observations revealed that the optical 
afterglow had a major rebrightening starting at ~ 0.01 days and ending at ~ 0.03 
days, which was followed by an initial power-law decay. Here we show that this 
early emission feature can be interpreted as the reverse shock emission superposed 
by the forward shock emission in an interstellar medium environment. By fitting 
Q\ | the observed data, we further constrain some parameters of the standard fireball- 

shock model: the initial Lorentz factor of the ejecta 70 > 120, the magnetic energy 
fraction e# > 4 x 10 -6 , and the medium density n < 2 cm -3 . These limits are 
I consistent with those from the other very-early optical afterglows observed so far. 

In principle, a wind environment for GRB 050525A is disfavored. 
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Introduction 



The gamma-ray burst (GRB) 050525A is a bright brief flash of gamma-ray radiation 
detected by the Swift Burst Alert Telescope (BAT) on 25 May 2005 00:02:53 (UT) (Band et 
al. 2005). It showed two peaks, with a duration of T90 = 8.8 ± 0.5 s in the 15-350 keV band 
(Markwardt et al. 2005; Cummings et al. 2005). Fitting to the Band spectral model yields a 
low-energy photon index ofa = 1.0±0.1 and a peak energy of E p = 79 ± 4 keV (Cummings 
et al. 2005). The fluence in the 15-350 keV band is (2.0 ± 0.1) x 10~ 5 ergs cm' 2 . This burst 
was also detected by other on-board instruments such as INTEGRAL and Konus-Wind, with 
peak fluxes of ~ 3.2 x 10~ 6 ergs cm" 2 s _1 (Gotz et al. 2005) and ~ 8.7 x 10~ 6 ergs cm -2 s _1 
(Golenetskii et al. 2005), respectively. The fluences were > 1.2 x 10~ 5 ergs cm -2 in 20-200 
keV (for the 12 second integration time) by INTEGRAL and ~ 7.8 x 10 _5 ergs cm -2 by 
Konus-Wind in 20-1000 keV (for the 11.5 second integration time), respectively. The time- 
integrated spectrum shows a peak energy E p = 84.1 ± 1.7keV by Konus-Wind (Golenetskii 
et al. 2005). About six minutes later, the ROTSE III telescope in Namibia was able to 
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obtain images, and detected an optical afterglow of magnitude 14.7 (Rykoff, Yost, & Swan 
2005). A presumed host galaxy was measured with a redshift of z=0.606, based on both 
[O III] 5007 and Hp emission and Ca H&K and Ca I 4228 absorption (Foley et al. 2005). 
Assuming cosmological parameters Qm — 0.27, Q\ = 0.73, and H = 71kms _1 Mpc -1 , the 
isotropic-equivalent gamma-ray energy is about 1.2 x 10 53 ergs. 

The follow-up observations of GRB 050525A revealed that the optical afterglow decayed 
as oc t~ L3 until about 0.01 days after the burst. Subsequently the afterglow showed a major 
rebrightening starting at ~ 0.01 days and ending at ~ 0.03 days, and then decayed as oc t~ L0 . 
This early-time emission feature is similar to the behavior of external reverse shock emission 
predicted by the standard fireball-shock model (Meszaros & Rees 1997, 1999; Sari & Piran 
1999a), which was first confirmed by the detections of early prompt optical emission from 
GRB 990123 (Akerlof et al. 1999; Sari & Piran 1999b). The reverse shock model was also 
suggested to interpret early optical emission of GRB 021004 (Fox et al. 2003a; Kobayashi 
& Zhang 2003a) and GRB 021211 (Fox et al. 2003b; Li et al. 2003; Wei 2003; Panaitescu 
& Kumar 2004b). Besides, the prompt optical-infrared emission from GRB 041219A was 
recently suggested to be the internal shock emission, since a strong correlation between 7- 
ray and optical-infrared signals was detected (Vestrand et al. 2005; Blake et al. 2005), which 
further confirmed the standard internal-external-shock model of GRBs. Mirabal et al. (2005) 
reported a break at ~ 0.4 days in the optical afterglow light curve of GRB 050525 A. In this 
paper, we show that the optical afterglow in ~ 0.4 days after the Swift trigger can be well 
understood due to the reverse shock emission superposed by the forward shock emission. 



2. Optical Emission from Forward and Reverse Shocks 

The standard model for GRBs and their afterglows is the fireball-shock model (for recent 
review articles, see Meszaros 2002; Zhang & Meszaros 2004; Piran 2005). In this model, 
the prompt emission of GRBs is ascribed to internal shocks and the long-term afterglow 
to external shocks. The prompt or very early optical emission is therefore explained as a 
consequence of the reverse component of the external shocks (Meszaros & Rees 1997; Sari & 
Piran 1999b). 

For shock-accelerated slow-cooling electrons with a power-law distribution, the typical 
spectrum of synchrotron emission is described by three power laws (Sari, Piran, & Narayan 
1998): (1) F v = {v/v m ) 1/z F^ max for v < v m ; (2) F v = (u /v m )-^l 2 F v , max for v m < v < u c ; 
(3) F v = {iy c /^m)~ {p " 1)/2 {^/^c)~ p/2 F U)max for v > v c , where u m , v c and F^ max are the typical 
synchrotron frequency, cooling frequency and peak flux, respectively. The parameter p is 
the index of the electron energy distribution. Here we ignore the synchrotron self-absorption 
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because its corresponding frequency may be much less than the optical band of interest. 

Both a forward shock and a reverse shock emerge when an ultrarelativistic cold GRB 
ejecta with initial Lorentz factor of 70 sweeps up a stationary cold interstellar medium: the 
forward shock propagates into the interstellar medium and the reverse shock propagates 
back into the ejecta (Katz 1994; Sari & Piran 1995). The emission of the forward shock is 
characterized by (Sari, Piran, & Narayan 1998) 



- - --° 1 1^)" 1/ Xl5feii)" 1/2 (^)" 3/2 (o^)«-V' 2 H, (2, 




where, t m j is the critical times when the break frequency, v m j crosses the observed frequency 
vr, ts is the fraction of the shock energy goes into the magnetic field, e e is the fraction of the 
shock energy goes into the electrons, g m = (p — 0.67) (p — 2) 2 / (p — l) 2 , g c = (p — 0.46)e~ L16p , 
Qmax = P + 0.14 (Granot & Sari 2002), E iso is the isotropic-equivalent kinetic energy, n is 
the density of interstellar medium in units of 1 cm~ 3 , z is the redshift of the burst, Dl is the 
corresponding luminosity distance, and is the observer's time in units of 1 day. The above 
equations are valid for p > 2. Afterglows with the hard electron spectrum of 1 < p < 2 have 
been discussed by Dai & Cheng (2001). 

For reverse shocks, it is possible to get a simple analytic solution in two limiting cases: 
thin shell and thick shell (which are corresponding to Newtonian and relativistic reverse 
shock, respectively; Sari & Piran, 1995). Since the reverse-shocked gas is separated with the 
forward- shocked gas by a contact discontinuity, which keeps the equality of pressures and 
velocities in both shocks, we can find the emission properties of reverse shock in the aid of 
the correlations between forward and reverse shocks at the crossing time t x (Kobayashi 2000; 
Kobayashi & Zhang 2003a). First, since the reverse- and forward-shocked gases have the 
same magnetic field and Lorentz factor at t x , the cooling frequency of the reverse-shocked 
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electrons is equal to that of the forward-shocked electrons, is c ,r{t x ) — v c ,f{t x ). Second, the 
electron's random Lorentz factor of the reverse shock is smaller than that of the forward shock 
by a factor of 7o/7 x (where 7 X is the Lorentz factor of both shocked regions at t x ), the typical 
synchrotron frequency of the reverse shock is by a factor of 7q/7 x l ess than that of the forward 
shock, v m ,r{tx) = (7o/7x) z/ m,/(^x)- This is basically valid for the thick shell case. As for a 
thin shell, it is more realistic with the expression, v m , r {t x ) = (734, x — 1) 2 /(7 X — l) 2 ^ m ,/(tx), 
where 734, x — (7 x /7o + 7o/7x)/2 is the Lorentz factor of the shocked ejecta in the rest frame 
comoving with the unshocked ejecta at the crossing time and 7 X ~ 7 /2 (Zhang, Kobayashi, 
& Meszaros 2003; Fan & Wei 2005). Third, the peak flux F v ^ max of an emission region is 
proportional to the electron number, the magnetic field, and the Lorentz boost. Since at 
t x , the electron number of the ejected shell is larger than that of the swept-up medium by 
a factor of 7 x /7o, we obtain F^ max , r (t x ) = ('~fl/lo)F V;max j(t x ). Therefore, we can directly 
write down the corresponding characteristic variables of the reverse shock emission, at t — t x , 




for the thin shell (i.e. Newtonian) case, and 
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for the thick shell (i.e. relativistic) case. Here, the GRB duration T is given by the shell 
width Aq/c, according to the internal shock model. As derived by Sari & Piran (1999a) 
and Kobayashi (2000), the typical synchrotron frequency and peak flux evolve as v m r oc t 6 , 
Fu,max,r oc t 3 ^ 2 at t < t x , and v m;r oc t -54 / 35 , F v>max ^ r oc £ -34 / 35 at t > t x , for the thin-shell 
case; v m r =constant, F v<max ^ oc t 1 / 2 at t < T, and v m r oc t~ 73 / 48 , F v ^ max ^ r oc t~ 47 ^ 8 at t > T, 
for the thick-shell case. Before the crossing time, the cooling frequency evolve as v c r oc t~ 2 
for the thin-shell case and u c , r oc t~ l for the thick-shell case; after the crossing time, u CyT 
turns into the cutoff frequency due to synchrotron cooling of electrons without acceleration, 
and it evolves the same as v m r for both cases. According to the emission properties shown 
above, we derive light curves of the emission from the forward- and reverse-shocked regions 
and thus can constrain the model parameters by fitting the early afterglow of GRB 050525A. 



3. Constraints on Parameters of GRB 050525A 

From the emission features described by equations (1), (2), (5), (6), (8), and (9), we see 
that the optical band is commonly between the typical frequency u m and cooling frequency 
v c in both forward and reverse shocks around the crossing time. Thus, the optical light 
curve of GRB 050525 A can be fitted by the superposition of two broken power-laws: (1) 
forward shock: F v j oc t l J 2 for t d < t m j, and F v j oc £ 3(1_p )/ 4 for t rf > t m j; (2) reverse shock: 
F^r oc i 3 ( 2p_1 )/ 2 for td < t x , and F v ^ r oc t^ d 27p+7 )/ 35 f or t d > £ x _ These scaling laws are 
valid for the thin-shell case, in which the reverse shock decelerates the shell insignificantly 
(Kobayashi 2000). 

In Figure 1 we present the R-band afterglow light curve of GRB 050525A between May 
25.0062 UT (At w 0.004 days) and May 25.5333 UT (At w 0.531 days). We find that, 
t m j 0.027 days, F u ^ max j ps 0.719 mJy, and p ps 2.3 can fit the light curve well for the 
component of forward shock emission. Given redshift z = 0.606 (Foley et al. 2005), the 
corresponding luminosity distance is D L = 1.1 x 10 28 cm. Using equations (3) and (4), we 
obtain 

E N 1/4 

e e ~8.3n^l- 

and 



e e ~ S.Sn 1 / 4 -£22- , (11) 
\10 53 ergs/ 



e B ~ 1.0 x lO- 5 ^ 1 ( — ) . (12) 
( 10 53 ergs/ 
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Furthermore, from the flux inferred by the light curve of the reverse-shock emission (depicted 
as a dashed line in Figure 1, where we see a flux F v>r = 0.75 mJy at t — 0.01 days), using 
equations (5) and (7), we further obtain 

(rp \ 1.658 x 0.825 / f \1.3 / ^ n n -2.965 

ife) (£0 (iFt) " a ' 67 (l^) (13 » 

By simple algebra from equations (11), (12) and (13), we obtain a constraint on the initial 
Lorentz factor of GRB ejecta, 




which is weakly dependent on the isotropic-equivalent kinetic energy E iso and interstellar 
medium density n. Assuming the efficiency of energy conversion is i] = 0.5, E iso is approxi- 
mately equal to the detected isotropic-equivalent 7-ray energy, i.e. 1.2 x 10 53 ergs. Thus, we 
can obtain a mutual constraint between e# and e e from equations (11) and (12). By elimi- 
nating the medium density n, we have e# = 4 x 10 _6 e~ 4 . Moreover, as a natural consequence 
of e e < 1, we have 70 > 116, €3 > 4 x 10~ 6 , and n < 2 cm -3 . These constraints are more 
clear than the cases of GRB 990123 (Sari & Piran 1999b), GRB 021004 (Kobayashi & Zhang 
2003a) and GRB 021211 (Wei 2003). Our preliminary fitting to the afterglow observations 
of GRB 050525A update to May 25.5333 UT (At « 0.531 days) is shown in Figure I. 

Additionally, we also give a constraint on the parameters for the thick-shell case. We 
can find from a fitting similar to that shown above that, the crossing time of the reverse 
shock, i.e., the GRB duration T is approximated by 

/ E iso \ °'° 6 f n \ ~ - 06 / 7 \ o- 2 5 
T ^ 553 V 10^3 ergs J ll^=sJ llO^J B ' (15) 

which is inconsistent with the observed T 90 ^ 10 s (Markwardt et al. 2005; Gotz et al. 2005). 
This implies that GRB 050525A seems to be a thin-shell case. 



4. Discussion and Conclusions 

The optical afterglow of GRB 050525A rebrightened starting at ~ 0.01 days and ending 
at ~ 0.03 days, which was followed by an initial power-law decay. In this paper, we have 
shown that this early emission feature can be interpreted as the reverse shock emission 
superposed by the forward shock emission in an interstellar medium environment. A good 
fitting to the observed light curve is shown in Figure 1. We further find the initial Lorentz 
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factor of the ejecta 70 > 120. Some other model parameters are also constrained: the 
magnetic energy fraction e# > 4 x 10~ 6 and the medium density n < 2 cm -3 . These limits 
are consistent with those from the other very-early optical afterglows observed so far (Sari 
k Piran 1999b; Kobayashi k Zhang 2003a; Wei 2003; Panaitescu k Kumar 2004b; Nakar 
k Piran 2005) and with the values of the forward-shock parameters (Panaitescu k Kumar 
2001). 

Our model is simplified under some assumptions: First, the reverse- and forward-shocked 
regions have the same electron and magnetic energy fractions, requiring that the ejecta is 
not initially magnetized. This requirement seems to be relaxed for GRB 990123 (Fan et 
al. 2002; Zhang, Kobayashi, k Meszaros 2003). The interaction of a magnetized ejecta 
with its surrounding medium has been investigated in details (Fan, Wei, k Wang 2004; 
Zhang k Kobayashi 2005). Second, we neglect the effect of inverse Compton scattering 
(ICS) on the reverse shock emission in §2. As shown by Wang, Dai, k Lu (2001a, b), this 
effect not only decreases the cooling frequency but also produces higher-energy emission. 
However, this decrease does not affect the limits derived above because the cooling frequency 
is much larger than the optical band at the times of our interest. Third, the surrounding 
environment is assumed to be a uniform interstellar medium, in which case two types of light 
curves are expected for the reverse-forward-shock emission combination: a rebrightening light 
curve (type I) and a flattening light curve (type II) (Zhang, Kobayashi, k Meszaros 2003). 
Early afterglows in wind environments have been studied by Chevalier, k Li (2000), Wu 
et al. (2003), Kobayashi k Zhang (2003b), Panaitescu k Kumar (2004b), and Zou et al. 
(2005). Generally, a flattening light curve in optical band is expected in the wind model 
(Kobayashi k Zhang 2003b; Zou, Wu, k Dai 2005), i.e., there is only type II light curve for 
the reverse-forward-shock emission combination. Since the early optical afterglow of GRB 
050525A shows a rebrightening (type I) light curve, the reverse-forward shock model in a 
wind environment for the early afterglow of GRB 050525A is not favored. 

The observed break at tj ~ 0.4 days in the optical afterglow light curve of GRB 
050525A (Mirabal, Bonfleld, k Schawinski 2005) may be due to an ultrarelativistic jet. 
This conclusion seems to be strengthened by Lulin Observatory around 17.7 hours after the 
burst (at which the R-band magnitude is about 21.3 ± 0.1, Chiang et al. 2005). If so, the 
jet's half opening angle 9 = 0.051(^/0.4 d) 3 / 8 (E iso /10 53 ergs)- 1 / 8 (n/0.1 cm- 3 ) 1 / 8 ( 7 ji r /0.5) 1 / 8 
(Sari, Piran, k Halpern 1999) and thus the beaming-corrected gamma-ray energy Ej et = 
1.3 x 10 50 (t j /0.4d) 3 / 4 (E iso /10 53 ergs) 3 / 4 (n/0.1cm" 3 ) 1 / 4 (77 7 /0.5) 1 / 4 ergs, which is by a factor 
of a few smaller than the mean energy release found by Frail et al. (2001) for n ~ 0.1 cm" 3 
and ?7 7 ~ 0.5. Even so, this energy still satisfies the relation of Ghirlanda et al. (2004), 
(£ ie t/10 50 ergs) = (1.12 ± 0.12) [(1 + z)E p /l00 keV] L50±a08 (see Dai, Liang, k Xu 2004). 
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Finally, what we want to point out is that besides the reverse-shock model discussed 
in this paper, the other plausible explanations of a rebrightening light curve include the 
density-jump medium (Dai & Lu 2002; Lazzati et al. 2002), pure Poynting-flux injection 
(Dai & Lu 1998), baryon-dominated injection (Rees & Meszaros 1998; Granot, Nakar & Piran 
2003; Nakar, Piran & Granot 2003; Bjornsson et al. 2004), two-component jet (Berger et 
al. 2003; Huang et al. 2003), neutron-fed fireball (Beloborodov 2003), temporal fluctuation 
of e e and/or e# and departure of the electron distribution from a power law (Panaitescu 
& Kumar 2004). Multiwavelength observations in the Swift era are expected to distinguish 
among these possibilities. 
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Fig. 1. — R-band afterglow light curve of GRB050525A. The fitting curve (solid line) is 
generated by superposition of forward shock emission (dotted line) and reverse shock emission 
(dashed line) with e B = 0.0025, e e = 0.2, n = 0.003 cm" 3 , E iso = 10 53 ergs, p = 2.3, 
7o = 200. These values are basically consistent with the thin-shell case. Data are taken from 
the GRB Coordinates Network (GCN): Malesani et al. (2005); Klotz, Boer, & Atteia (2005); 
Torii & BenDaniel (2005); Homewood et al. (2005); Rykoff et al. (2005); Yanagisawa, Toda 
& Kawai (2005); Milne, Williams, & Park (2005); Mirabal, Bonfield, & Schawinski (2005); 
Cobb & Bailyn (2005); Klotz et al. (2005). 



